Instanton modifications of the bound state singularity in the Schwinger Model 
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We consider the quark-antiquark Green's function in the Schwinger Model with instanton con- 
tributions taken into account. Thanks to the fact that this function may analytically be found, 
we draw out singular terms, which arise due to the formation of the bound state in the theory — 
the massive Schwinger boson. The principal term has a pole character. The residue in this pole 
contains contributions from various instanton sectors: 0, ±1, ±2. It is shown, that the nonzero ones 
change the factorizability property. The formula for the residue is compared to the Bethe-Salpeter 
wave function found as a field amplitude. Next, it is demonstrated, that apart from polar part, 
there appears in the Green's function also the weak branch point singularity of the logarithmic 
and dilogarithmic nature. These results are not in variance with the universally adopted 5*-matrix 
factorization. 

PACS numbers: ll.10.Kk, ll.10.St, 11.55.-m 



I. INTRODUCTION 

The Schwinger Model (SM) Q, describing a massless 
fcrmion field in interaction with a gauge boson in two 
space-time dimensions, has become a very fruitful sys- 
tem in field theory. Above all it possesses many of the 
features, one is convinced to find in the theory of strong 
interactions. Therefore, thanks to its relative simplic- 
ity, SM is expected to be a perfect testing layout for 
the various nonperturbative aspects of QCD. They in- 
volve confinement, topological sectors, instantons, con- 
densates. Other interesting and nontrivial features one 
should mention, are the existence of anomaly and gauge 
boson nonzero mass generation without the need of intro- 
ducing an auxiliary Higgs field. Naturally some of these 
properties are related with one another being in fact dif- 
ferent manifestations of the same physical phenomena. 

Another aspect, that should be appreciated in SM, 
is the exact solvability of the equations, which allows 
one to perform many calculations without unpleasant 
approximations and additional assumptions or simplifi- 
cations. One example is the quark propagator (the term 
'quark' means here the fundamental fermion, due to its 
confinement property) , which has already been found by 
Schwinger, but also other functions stand wide open. For 
instance, in our previous work |2|, we were able to ob- 
tain the explicit formulae for the whole set of four-point 
Green's functions. 

A very troublesome and long-standing question in 
Quantum Field Theory is the consistent description of 
bound states. The problem becomes particularly severe 
when such a state has a relativistic character and re- 
tardation effects in the interactions among constituent 
particles cannot be neglected. Over half a century 
ago, the equation for the bound state wave function 
- the so called Bethe-Salpeter equation — was formu- 
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lated 0,0). This equation is, however, extremely com- 
plicated. Firstly, it is a multidimensional integral equa- 
tion and these are always more challenging than differen- 
tial ones. Secondly, being an equation for the unknown 
bound state wave function, it requires, as an input, the 
knowledge of the nonperturbative propagators for ingre- 
dient particles and the interaction kernel between them. 
Naturally no such quantities are known in realistic field 
theories. From the very beginning one is then doomed to 
extremely strong approximations. As to the calculations 
in field theoretical models, there are only few examples 
(for instance the so-called Wick-Cutkosky Model [5|), 
where it is possible to find solutions but even in that case 
one is forced to simplify the equation either by neglecting 
relative time dependence or simplifying the interaction 
kernel together with one-particle propagators. 

The other possibility in investigating bound states, in- 
stead of solving Bethe-Salpeter equation, is to consider 
four-point Green's function in the appropriate channel, 
and find the residue of the pole corresponding to a given 
bound particle @, Q. This residue is constructed just 
from the Bethe-Salpeter wave functions. However, know- 
ing the full four-point function is usually as unlikely as 
solving exactly the Bethe-Salpeter equation itself. One 
should then all the more appreciate the SM, for which 
the exact form of the Green's function in question, was 
found 0, [|| . This is the main reason, why we decided 
to look for the bound state singularity in this model, to 
analyze the nature of this singularity, and to derive the 
formula for the Bethe-Salpeter function. The calculation 
performed primarily in section IIIII shows the complete 
agreement with general considerations [fj, [7|, [9( . 

A certain new aspect appears however, if one takes into 
account that the SM is a topologically nontrivial the- 
ory 0, El E El El El - As it is well known, fcrmion 
Green's functions gain contributions from nonzero in- 
stanton sectors. This is connected with the existence 
in the theory of the (infinite) set of topological vacua, 
the superposition of which constitutes the true ground 
state. The expectation values for field operators in this 
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new state, become now sums of matrix elements taken be- 
tween various topological vacua. If a given operator has 
non-vanishing off-diagonal elements, then they have to 
enter into the calculation, apart from the main, diagonal 
contributions. Distinct topological vacua correspond to 
different, and topologically inequivalent, configurations 
of the gauge field, and the corresponding transition be- 
tween them is described by the instanton field. There 
is a proportionality between topological index of a vac- 
uum and the chiral charge and, therefore, off-diagonal 
elements may appear only for the chirality variant oper- 
ators. To that category belong the products of fermion 
fields and, consequently, we may expect the appearance 
of certain new terms in the Green's functions involving 
quarks. The explicit calculations for the two- and four- 
point functions has been done in [l6l |. where these sup- 
plementary terms were explicitly given. 

Now the question arises if, and how, these contribu- 
tions modify the behaviour of the two-fcrmion Green's 
function close to the bound state singularity. To this 
question we devote section TlV Al (one instanton contribu- 
tions) and section [IV Bl (two instantons contributions). 
Our results indicate two modifications. Firstly, we see 
that the residue in the polar term acquires new ingredi- 
ents: the Bethe-Salpeter wave function is modified — and 
this is quite obvious, since it is defined through the prod- 
uct of the two quark fields, which is sensitive to instanton 
background — and the other term appears, which makes 
the residue nonfactorizable. Yet, it may still be written 
as a sum of several factorizable terms. Secondly, apart 
from the pole, we find terms, which show another type of 
singularity — branch point singularity of the logarithmic 
and dilogarithmic nature. These contributions are small 
in comparison with the principal term, and in the usual 
bound state pole derivations d, 0] , are neglected by the 
proof construction itself. In the SM, we happily dispose 
the full Green's function without being obliged to make 
any addition assumptions or simplifications, which usu- 
ally remain beyond our control, and those terms survive. 

On the other hand, the calculated instanton con- 
tributions disappear for transition amplitudes between 
asymptotic states and to that extent our results are 
in agreement with the analyticity properties of the S- 
matrix [|, [U H El US [SH, 112. For the confined parti- 
cles one cannot require the cluster decomposition on the 
level of the Green's function. The main argument refer- 
ring to the possible separation of certain subgroups of 
particles into the distant space-time regions, may be ap- 
plied only for hadrons, but neither for individual quarks 
nor even qq structures, if they do not bear the complete 
set of quantum numbers of a meson. It is then not sur- 
prising, that the quark Green's function reveals much 
richer analytical structure, than the hadronic S'-matrix 
does. 

Finally, in section [V] we calculate directly the Bcthc- 
Salpeter wave function as a qq field amplitude between 
one-particle state of momentum P and the vacuum state, 
and find that it corresponds to that obtained from the 



residue of the factorizable part of the preceding section. 

II. BASIC DEFINITIONS AND NOTATION 

The Lagrangian of the Schwinger Model, describing a 
fermion (quark) field '3/ in interaction with a gauge boson 
A 11 , has the following form: 

C(x) = *(x) (i^d„ - g^A^x)) *(x) 

- \F^{x)F llv {x)-^{d^{x)f , (1) 

where A fixes the gauge and g is the coupling constant. 
The space-time has one time and one space dimension. 
In this world, the most simply is to choose the quark 
field in the two-component form. In this case, the 2x2 
Dirac gamma matrices may have the representation given 
below: 

For the metric tensor we choose the convention: 

.? 00 = -g 11 = i , 

and the antisymmetric symbol is defined by: 

e 01 = -e 10 = 1 , e 00 - e 11 = . (3) 

All the Green's functions may be obtained by differenti- 
ating over external currents of the generating functional, 
which is, as usually, defined by the Feynman path inte- 
gral: 

Z[ V ,% J]= J V^V^VAe l I d2x ^ + ^ +J ' lA ^ , (4) 

As we have already written in the Introduction, this the- 
ory is topologically nontrivial. The true vacuum, which 
is usually called the ^-vacuum, becomes a superposition 
of the infinite set of degenerate vacua, bearing different 
topological indices n [H HH : 

oo 

\9>= em9 \ n > ■ ( 5 ) 

n— — oo 

In this new vacuum, the generating functional ([4]), which 
constitutes the Feynman form of < 6\6 > v ,rj,j, breaks 
into pieces: 

oo 

Z[ V ,rj,J} = £ e lk9 Z^[ V ,%J}. (6) 

k— — oo 

In this formula — contrary to © — the summation runs 
over instanton number, i.e. k corresponds to the change 
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of the topological indices in the two vacua in the expres- 
sion: < 6\6 >„m,j- Henceforth we will use the com- 
mon symbol |0 >, for denoting the true vacuum. In each 
of Z( fe )'s the functional integration with respect to the 
gauge field , is restricted to a specific class of configu- 
rations bearing the Pontryagin index k. 

Now the qq Green's function may be represented as the 
appropriate functional derivative of Z: 

G a b,cd(xi,X2;X3,Xi) = (7) 
=< 0|T(* a (xi)* 6 (x 2 )* c (a; 3 )* d (x4))|0 >= 
6 4 Z 



= N 



5i] a (xi)5r) b (x2)5r) c (x 3 )Sr) d (x i ) 



where N is a constant ensuring that < 1 >ri.rj.j=o= 1- 
In the following, we are only interested in singularities 
that may be associated with bound states, and hence, 
throughout this work, by G we understand only the con- 
nected part of this function. 

Due to the fact, that Z expands in the sum of 
the same refers also to the Green's function. However, 



this sum is now limited only to the low values of k. Con- 
trary to the expression exp {i J d 2 x \jp$ + ^frj\ }, which 
contains arbitrary high powers of \& and ^ and, there- 
fore, contributes to all sectors in ([6]), the four-fermion 
function gets terms with at most k — ±2. This is due to 
the fact, that four fermion fields can change the chirality 
of the state by (at most) 4, and the topological vacuum 
index is a half of the chiral charge [23j . In the following 
sections we investigate in detail all possible terms. 



III. QUARK-ANTIQUARK BOUND STATE 
WITHOUT TOPOLOGICAL CONTRIBUTIONS 

In this section we will concentrate on the O-instanton 
term of the function G a b,cd{x\, x 2 ; x 3 , £4). It is a great 
merit of the Schwinger Model, that this function may 
be found exactly, in an explicit and analytic way. The 
details of our calculation have already been given in 
and we will not elaborate on them here, but instead recall 
only the final result: 



G l°bU {xi,x 2 ;x 3l x 4 ) = 



= 2^ Sac ^ Xl ~ x z)Sbd{x2 - Xi) + [S(xi - x 3 )j 5 ] ac [S(x 2 - x 4 )7 5 ]fcd} e-xp{ig 2 [(3{xi - x 2 ) - (3{xi - x t ) 
-f3(x 2 - x 3 ) + /3(x 3 - x 4 )}} + ^{S ac (x 1 - x 3 )S bd (x 2 - x 4 ) - [S(xi - x 3 )Y J ]ac[S(x 2 - x 4 )7 5 ] h(i } 



x exp{-igi 2 [/3(a;i - x 2 ) - /3(xi - £4) - f3(x 2 - x 3 ) + fi(x 3 - £4))]} - 

I 



£3 
c 



Xi 

d 



(8) 



where the index (0) refers to the instanton number. The 
function (3(z) is defined through the two-dimensional mo- 
mentum integral: 



d 2 p 



1 - e lpz 



(27r) 2 (p 2 - g 2 /-K + ie)(p 2 + ie) 



(9) 



and may be expressed in terms of MacDonald or Hankcl 
function, depending on whether z is space- or time-like, 
but the above representation is more suitable for us. We 
recall that in 2D the coupling constant g is dimensionfull 
and the quantity /x 2 = g 2 /ir plays a role of the mass 
squared of the so called Schwinger boson - a qq bound 
state we are interested in. 

The object S(x) in © is the full quark propagator, 
which has the known form: 



S(x) = <S (a:) ex P [- i 9 2 P( x )] 



(10) 



with Sq(x) 



2tt x 1 



being the free propagator. The 



term 'full' used above, means here 'full in the O-instanton 
sector', i.e. it is given in the form of the original 
Schwinger works [l[. However S - bilinear in 'S's - 



must later get contributions from the nontrivial topolog- 
ical configurations 16]. 

Now we will concentrate on the structure of the qq 
function (|8|) and try to reveal, whether it has a pole in 
the t-channel. To this goal we introduce new space-time 
variables: X, Y, x and y, defined in the following way: 



X = ~(xx + x 3 ) , x = xi-x 3 , 



Y = ~(x 2 +x A ) , 



y = x 2 - x 4 



(11) 



Thanks to the translational invariance, which manifests 



itself through the dependence of G on Xi 



's only, we 



have three independent variables: Z = Y — X, x and 
y. From the four terms present in l[5]). the two explic- 
itly written are of interest for us, since the two other, 
obtained by the antisymmetrization, contribute to the 
singularities not in t, but eventually in the u-channel. 
First of the two interesting terms, denoted by /, takes 
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the form: 



G^\Z;x,y)=\{S{x) ® S(y) + [5(x) 7 5 ] ® [S(y) 7 5 ]} 

x exp{i 5 2 [/3(-Z + (x - y)/2) - + (x + y)/2)) 

-0(-Z -(x + y)/2)) + (3(-Z -{x- y)/2))]} , (12) 

where we made use of the evenness of f3: (3(x) — f3(—x). 
For brevity we omitted the spinor indices, and used the 
tensor product notation, which should not cause any con- 
fusion. The pole corresponding to the bound state should 
be found in the complex plane of P 2 , where P is the two- 
momentum canonically conjugated to Z . The matrical 
part depends only on the relative coordinates x and y and 
enters into the residue of this pole, but does not influence 
its position. Let us then perform the Fourier transform 
over Z of the exponent containing (3 functions. Making 
use of © we my write the quantity in question as: 



iPZ -4*9* J ^e- ipZ h(p,x,y) 



(13) 

where h(p,x,y) = l^^i+^^l^ ■ Expanding the ex- 
ponent in series, we can execute the Z integration ob- 
taining: 

F\P,x,y) = £ t^pl In{P:X>y) + {2 n) 2 S^(P) , 

* — ' TV. 

71=1 

(14) 

The I n 's satisfy the following recurrent relation: 

In (q,x,y) = (15) 
f d 2 p sm(px/2)sm(py/2) 
J (2i:y (p z - p, 2 + it){p z + it) 

It is obvious that Ii(p,x,y) has a pole at p 2 = /i 2 . The 
other singularity is at p 2 = 0. But what about terms of 
higher nl Do they alter the behaviour at p 2 = /i 2 ? Let 
us first construct ^(p, x,y). According to (fT5|) it has the 
form: 



h (q,x,y) 



(16) 



cPp sm(px/2)sin(py/2) 
(2tt) 2 (p 2 - [i 2 + it)(p 2 + it) 



h(q-p,x,y) 



The possible singularities in q 2 , exhibited in 1%, can come 
from the singularities of the integrand functions. There 
are four of them: obliviously p 2 — 0, p 2 = /1 2 , and the 
other two originating from I±. These are: (p—q) 2 = and 
(p — q) 2 = fi 2 . It is not possible to satisfy all these four 
or even three conditions simultaneously, so the singular 
points may appear only as a result of the merging of two 
singularities. In accordance with Landau procedure @, 
HEHil , we introduce two parameters (a and j3) satisfying 
a, (3 > and a 2 + (3 2 ^ 0. With the assumption ap + 
(3{p — q) = 0, the following possibilities occur: 

1. ap 2 — and (3{q—p) 2 = 0. If one of the parameters 
equals zero, then we get p = 0orp — q = and the 



eventual singularity is cancelled by the sine func- 
tions present in the numerator in (|16p . either by 
that explicitly written or by the one hidden in the 
formula for I\. If both parameters are nonzero, on 

can easily obtain q 2 = (^-jf-^J P 2 = 0. Thus, there 
may only arise a contribution to the singularity at 
q 2 = 0. It cannot, however, be worse than that in 

2. a(p 2 — (i 2 ) = and (3{p—q) 2 =0. If a — 0, we have 
q = p and, in general, one might get a contribution 
to the pole at q 2 = fi 2 (actually it would become 
a branch point), but it does not happen because of 
the sine functions becoming zero for q — p = 0. If 
(3 = 0, we obtain p = 0, which stays in contradic- 
tion with the assumption p 2 = /1 2 . For a, f3 7^ 0, 
we find a 2 p 2 = f3 2 (p — q) 2 , which is incompatible 
with p 2 = [i 2 and (p — q) 2 = 0. 

3. ap 2 — and (3((p — q) 2 — /i 2 ) = 0. This case is 
almost identical with the previous one and need 
not be separately considered. 

4. a(p 2 ~fi 2 ) — and (3[(p- q) 2 - [i 2 ] =0. If one of the 
parameters is zero, then we find a contradiction of 
ap + (3(p— q) = 0, with either p 2 — /i 2 or (p — q) 2 = 
[x 2 . What remains is to assume a, j3 7^ 0. If so, 
we have a 2 p 2 = f3 2 (p— q) 2 and after cancelling the 
arising factors fi 2 on both sides, we get a = (3. This 
leads to q = 2p and a singularity emerges now at 
q 2 = 4/z 2 . This is the two-particle singularity. 

Summarizing, we observe that in 7 2 the pole at [i 2 disap- 
pears and a new singularity arises at (2/i) 2 . Now we are 
in a position to carry out a proof by induction: the as- 
sumed singularities under the integral in (|15[) are: p 2 = 0, 
p 2 = fi 2 , {p-q) 2 = and {p — q) 2 — ((n— l)^) 2 . We have 
again four possibilities, the analysis of which is identical 
as above, except the last point, where we now get the 
equation a 2 — (n — 1) 2 (3 2 . This leads to q = np and 
q 2 = n 2 ^ 2 . Hence, the conclusion is, that the singularity 
corresponding to n intermediate particles (q 2 — n 2 p, 2 ) is 
present only in /„. This means, that the pole at fx 2 may 
be simply picked up from the first term of the expan- 
sion (fit]) . 

Now we need to apply the same to the second expres- 
sion in the Green's function ([5]). Again we take out the 
exponent and perform the Fourier transform, the result- 
ing expression being called F n (P,x,y). Happily both 
exponents differ only by the sign, so we are able to im- 
mediately write down the final result. The total contri- 
bution (i.e. I + II) to the bound state pole in G^ ' is: 



-4iir [5(x) 7 5 ] <g> [S(y)j 5 ] 
sm(Px/2)sm(Py/2) 



P 2 



(17) 



We see the perfect factorization of the residue, and the 
quantity: 

3>p{x) = -2y/^S(x)-f 5 sin(Px/2) , (18) 
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plays a role of the Bethe-Salpeter amplitude [1, H, 0] or 
the wave function. The minus sign is chosen for later 
convenience. Note, that this quantity is found explicitly 
and exactly, without any simplification referring neither 
to the relative time x° nor to the structure on the in- 
teraction kernel. It is an essential point, since even in 
the simplified model theories, one usually is limited to 
the case of the so called 'equal time wave function' (i.e. 
defined only on the hypersurface x° = in the center 
of mass frame), which neglects retardation effects and 
consequently leads to the nonrelativistic bound state, or 
ladder approximation for the interaction kernel. On the 
other hand, the function (fT8)) is a fully relativistic object. 



dom, corresponding to the transition amplitude in ques- 
tion, is then zero and the tunnelling phenomenon disap- 
pears. This does not mean, however, that the notion of 
(9-vacuum is useless, since topological vacua (contrary to 
the vacuum) violate cluster property (3(j. We are the 
fated to use ([5|) and must put up with modification in 
the fermion Green's functions 

The contributions to the four-point function were al- 
ready found in our previous work (l6| . Below we will 
concentrate on the eventual bound state pole in the for- 
mulae. 



IV. TOPOLOGICAL CONTRIBUTIONS TO 
THE QUARK-ANTIQUARK BOUND STATE 

As it has been mentioned, the vacuum in the Schwinger 
Model is nontrivial: it is a superposition ([5]) of vacua 
corresponding to various topological numbers. As a con- 
sequence, all chirality changing operators — and to that 
class the product of fermion fields belongs — which have 
nonvanishing matrix elements between different topolog- 
ical vacua, acquire contributions from nonzero instanton 
sectors. We recall here, that if the theory contains mass- 
less fermions (at least one) — and this is just the case 
of the Schwinger Model — the amplitude of the spon- 
taneous tunnelling between different topological vacua 
vanishes [23L |27| . Mathematically it is expressed through 
zero value of the euclidean Dirac operator determinant: 
det[i (jt— eA\, because of the appearance of zero eigenval- 
ues when bears the instanton index [24 1271. l28l. 129 . l3l| . 
The functional integral over fermionic degrees of free- 



A. 1-instanton sector 

By virtue on the Atiyah-Singer index theorem [3, H3, 
HH, 0, HH , the massless Dirac operator in the external 
field bearing nonzero instanton number k has zero modes. 
In the present section we consider the case of k = ±1, 
for which there is only one such mode. Its detailed form 
may be found in the literature [1(1 [H, HH) and we will 
not recall it here. This mode gives now its contribution 
while calculating the Feynman path integral Since 
it corresponds to the zero eigenvalue, it disappears from 
the quadratic part of the Lagrangian and remains only 
in source terms. Due to the Grassman character of the 
quark fields, of the four differentiations in ([7]), two must 
be saturated by sources multiplying zero modes. The 
other two are performed in a common way. The result 
for the four-point function, with both k = +1 and k = — 1 
instanton numbers taken into account, obtained in [ll| . 
is: 



G {1) (x 1 ,x 2 ;x 3l x 4 ) = -^J/je 7B (cos6» - ij 5 sin<9) <g> S (x 2 - x A ) [(1 <g> 1 - 7 5 <g> 7 5 ) 

xe ig 2 [P(xi-x a )-f3(x2-X4,)~f3(x 1 -x i )+f3{x 1 -X2)+f3(x2-x 3 )-f3(x 3 -x i )] |(l0l^.^ 5 g ,-y 5 ) i 
x ^a 2 mxi^x 3 )-f3(x 2 -x i )+f3(x 1 -x^-f3(x 1 -x 2 )^p(x 2 -x 3 )+f3(x 3 -x i )]^ + antisymmetr i zation . 

I 



The antisymmetrization runs here both over 'initial', 
i.e. x\ and x 2 and 'final' i.e. x 3 and X4 (accompanied 
by the appropriate spinor indices change) variables in 
G{x\, x 2 ; X3, X4). Thus p9[) gets three more terms, which 
are not explicitly written here for brevity. One of them 
(with the simultaneous substitution 1 <-> 2 and 3^4) 



again contributes to the pole in the i-channel, the other 
two only in the it-channel and will be omitted. Introduc- 
ing center of mass and relative variables (fTTj) . similarly 
as it was done in section Hill we obtain for the first term 
in CEU): 



G (1)/a (Z; x, y) = --^e~> E [(cos0 - i-f 5 smdy g2p{x) <g> S(y)] (1 <g> 1 - 7 5 <g> 7 5 ) (20) 
x exp{ig 2 [-p(-Z +{x + y)/2) + 0(-Z + (x - y)/2)) + 0(-Z - (x + y)/2)) - 0{-Z - (x - y)/2))]} , 
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The last exponent in the above, per analogiam to (|13p . 
may now be written in the form: 



dfp -i P Z 



(21) 



where Ji(p,x,y) = [°^l%^ 2 +i£) ■ Now we mi S nt in- 
troduce J„'s and consider the recurrent relation similar 
to (|15[) . but it is unnecessary. Due to the presence of 
the sine function in the numerator, we can directly apply 
the analysis of section IIIII and come to the result that the 
contribution to the bound state comes only from the first 
term of the series for the exponent (f2"Tj) . As to the second 
term (in our nomenclature it would be called lb) in (|19p 



the exponent has only the sign inverted and need not be 
elaborated on. Their joint effect is then simply to cancel 



1 <g> 1 in (HU) and to double 7 5 
to the bound state is therefore: 



7 5 . Their contribution 



cos(Px/2)sm(Py/2) 



P 2 - /i 2 



(22) 



where, instead of cos 8 — ij 5 sin 8, we have written e - * 7 6 . 
The antisymmetrized term, which is left so far, has an 
analogous form, but with x and y swopped, and tensor 
products inverted. Our final result in this sector is: 



x cos(Px/2) sin(Py/2) + [S(x)j 5 ] <g> [e^^e" 47 ^ 5 ] sin(Pa;/2) cos(Py/2)} . (23) 



Thanks to the quark field being massless, which results in 
the chiral invariance of the initial Lagrangian, the depen- 
dence on the parameter 8 may be gauged away, so one 
might simply put 8 — 0. Comparing the two obtained 
equations (|17[) and (|2"3"]) we see that introducing a new 
Bethe-Salpeter amplitude in the form: 



As we shall see in the following section, this is actually 
the case, but it contributes also other terms that modify 
the factorization and even change the character of the 
singularity at P 2 = /i 2 . 



-2V^S(xb 5 sin(Pa;/2) 



(24) 



^ -7 B ^9 2 /3(^-»7 5 7 5 cos ( P;E / 2 ) 



20? 



the residue factorization in the pole will be restored i.e. 
one will have: 



Gb.s.(P;x,y) 



$ P (s) ® $ P (y) 
P 2 - /J, 2 + ie 



(25) 



if the 2-instanton sector contributes to the bound state 
a term (with the appropriate constant coefficient): 



[e i9 ^e - * 7 e 7 5 ] 



[ e V/3(y) e -«7 5 e 7 5j 



x cos(Px/2) cos(Pj//2) 



2-instanton sector 



In the 2-instanton sector there are two zero modes of 
the the Dirac operator [13, EE E3, IH HI ■ Now all of the 
differentiations in ([7J are saturated with sources standing 
with zero modes. The 2-instanton contribution to the 
Green's function, which was evaluated in [l6j], has the 
form: 



4 



( Xl ,xr,x 3 , x 4 ) = -^le 47 - [e™ (l - 7 5 ) ® (l - 7 5 ) (4 + x\) (26) 

x (-x° 4 + x\) + e- 2W (1 + 7 5 ) ® (1 + 7 5 ) {-x\ + x\){xl + x\)} exp {ig 2 [p(xx - x 4 ) 

+(3(x2 — X3) + (3(xi — X2) + 0(x3 — X4) + (3(x\ — x%) + (3{x2 ~ 2:4)] } + antisymmetrization . 

Using again the variables Z, x and y, performing the required antisymmetrization, and replacing f3 functions with 
integral formula ©, we obtain, with the additional use of the identity e~ iei (1 ± 7 s ) = e Tl9 (l ± 7 s ): 
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G (2) (Z;x,y) 



a 4 
64tt 2 



J9 0iv)p—il 



2 ^(x-y) 2 h 5 ®7 5 



+ 1(Z 2 - 1(1 - y) 2 )((l - 7 5 ) ® (1 + 7 5 ) + (1 + 7 5 ) ® (1 - 7 5 )) 



(27) 



1 



a/3 Z Q (x - ® 7 5 + 7 5 ® 1)] exp [iwi/i 2 



d 2 p 1 



cos{px/2)cos(py/2)- 



(2ir) 2 (p 2 + ie)(p 2 - fi 2 + ie) 



with e Q/3 being the antisymmetric symbol defined in |3|). 
Now we will show that the first term (i.e. Z 2 times 
7 5 ® 7 5 ) gives just the contribution needed to restore 
the factorization l|25[) . 

One could suspect that there might be some cancel- 
lations, among various expressions in the formula ([27]) . 
Below we shall argue that the cancellation that would 
lead to complete disappearance of any term, is not pos- 
sible due to their independent matrical structure. Only 
cancellations among a specific matrix elements — but not 
matrices as a whole — are admissible. Consider namely 
the expression: 

a 7 5 ® 7 5 + 6[(l-7 5 )® (l + 7 5 ) + (l+7 5 )® (1"7 5 )] 
+ c(l <£> 7 5 + 7 5 ® 1) (28) 

and multiply it firstly by (1 + 7 5 ) <8> (1 + 7 5 ) and secondly 
by (a-7 5 )®(a-7 5 ). We obtain (a + 2c) (1 +7 5 ) ® (1 + 
7 5 ) + 6-0 and (a-2c)(l-7 5 )<8> (l-7 5 ) + 6-0 respectively. 
On the other hand if we take the double trace (i.e over 
first matrix and over second matrix in the tensor product 
indices), we get a-0 + 8- b + 0-c. This shows that the 
three tensors listed in (f2"8"]) are really independent. 

Let us observe now, that the quantity in the last expo- 
nent of (|2"7| , after rearranging the terms under the inte- 
gral (the integrand function will be called below M) and 
partially executing it, may be rewritten in the following 
way: 

4 ™^ 2 / i^M(Z,p,x,y) = -2 lE -\n(-LL 2 Z 2 /4) 



(2tt) 2 
d 2 p 



-ipZ 



i 



p 2 — \L 2 



(29) 



(2tt) 2 

d 2 p _ ivZ 1 - cos(pa;/2) cos(py/2) 
(27r) 26 (p 2 + ie)(p 2 -Li 2 + ie) ' 

If we expand, as before, the exponent of the last two 
terms and take the Fourier transform over Z, it be- 
comes obvious that both of them contribute to the pole 
at P 2 = /j, 2 . The further analysis is again identical to 
that of Sec. IIII1 except for the difference, that now there 
is additional cancellation between the integrals in ([29]) . 
Finally all the factors combine together in (f2"T)) and we 
obtain the expected term: 



47T 

cos(Px/2)cos(Py/2) 



This is exactly what we needed for ([23)) . This is not, 
however, the whole story. The similar contribution comes 
from the second matrical term in (|27[) . which also con- 
tains Z 2 . Since we have already saturated (|25[) this sup- 
plementary part violates the factorization. It has the 
form: 

■ 2 

Zlt c 2lEf c tg 2 f3(x) c -te~, 5 ^ e ig 2 P(y) e ~i9^ v 
87T 

x [(1 - 7 5 ) ® (1 + 7 5 ) + (a + 7 5 ) ® (a - 7 5 )] 

cos(Pir/2)cos(P?//2) 



P 2 - n 2 



it 



(31) 



(30) 



Due to the above mentioned lack of cancellations among 
independent matrix structures, this term must survive. 
Moreover it cannot cancel with the (x — y) 2 term with 
the same matrical coefficient since, contrary to (f3~Tj) . the 
latter disappears for x = y. Besides, it does not con- 
tribute to the bound state pole at all. This can be 
easily argued as follows. After rearrangements identi- 
cal to (|29p one gets the additional 1/Z 2 coming from the 
logarithm in the exponent. This can be neutralized, if 
we act with the two-dimensional d'Alambert operator in 
variable P. When this is done, we obtain a first order 
pole at P 2 = ii 2 . Now the question is: given a func- 
tion / satisfying: dpf(P) ~ pgjn ■ Does / also dis- 
play a pole? The answer is negative. It may easily be 
found that f(P) ~ Li2(P 2 /fJ. 2 ) and dilogarithm has a 
limit 7r 2 /6 when its argument goes to unity. Therefore 
it may be neglected in comparison with the polar term. 
One should emphasize, however, that although diloga- 
rithmic function has a well defined limit, it is not an an- 
alytic function in this point. It is a branch point and the 
function has a cut on the half line [/i 2 , oo[. The presence 
of cos(Px/2) cos(Py/2) does not change the conclusions, 
since cosine is an analytic function, and may only make 
things better. On the other hand for x = y = cosines 
can be replaced with unities and a branch point arises. 
Such terms are usually lost in the common approach, 
where one fishes out only the polar contribution of the 
bound state. 

For the term proportional to e a "Z a (x — y)p, similar ar- 
gumentation is applicable, with a difference that we now 
have to do with logarithmic function, which again may 
be neglected if compared to the polar term. Hence, as we 
see, due to the topological background, not only nonfac- 
torization appears, but the character of the singularity 
in P 2 = /j 2 is slightly changed. 
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These results do not, however, contradict the 5-matrix 
pole factorization [U (H [ij, Hfl, . The S-matrix de- 
scribes the transition amplitudes among physical asymp- 
totic states. In the Schwinger Model these are the Fock 
states of Schwinger bosons. Consequently we are limited 
only to the current-current Green's functions and this re- 
quire putting x, y — > 0, and taking the trace over spinor 
indices with 7^ ® 7^ matrices. Since the trace of an odd 
number of 7's is zero, it is easy to observe that the only 
contribution to the pole comes now from (|17[) and factor- 
ization reappears This is what should be expected 
since the vector current does not couple to instantons. 

Summarizing, that part of the full qq Green's function, 
in which one is limited to the polar terms, has the form: 

iG b . a XP;x,y)= 1 [$ P (x)®$ P (y) (32) 

Jr z — fi z + le 

+¥ p \x) <8 * { p\y) + ® t> { P \y)} , 

with 

9&>(x) = -V—e~< E e l92 ^ x) e- l9 ~<\t± 1 5 ) , (33) 
2v 2n 

and $p(x) given by (f24|) . 

V. WAVE FUNCTION AS A FIELD 
AMPLITUDE 

The question arises, how the obtained results harmo- 
nize with the Bethe-Salpeter wave function defined as an 
amplitude of fields. Fortunately in the Schwinger Model 
we have all tools, which allow us to find it explicitly even 
in the nontrivial topological sectors. Let us then consider 
the amplitude in question: 

Xp(xux 2 )=<0\T(y(x 1 )y(x 2 )\P> . (34) 

The calculation of this quantity turns out to be very sim- 
ple, if one already knows the expression for the full quark 
propagator in the instanton background [16| : 

Sfull(z) = S(x) + (cosfl - if sinfl) e T*+» 2 0(*) . 

(35) 

The state \P > is the one Schwinger boson state of two- 
momentum P, and (|34p is expressible, thanks to the for- 
mulae of the LSZ formalism, through the vertex function. 
It is known that, thanks to the two gauge symmetries 
(U(l) and Ua(X)), the vertex function may in turn en- 
tirely be constructed from the quark propagator. This is 
a miracle of the Schwinger Model. One then gets: 

Xp(x 1 ,x 2 ) = iy/ir[e~ lPxi j 5 S(xi - x 2 ) 

+e- iPx *S(x! - x 2 )j 5 ] . (36) 

If we now substitute for the quark propagator the ex- 
pression Sfuii, we obtain the formula for the amplitude 
Xp in all instanton sectors in question: 0, ±1 (recall that 



we now have to do with an object billinear in quark 
fields). We now pass to the CM and relative coordi- 
nates: x\ = X + x/2 and x 2 = X — x/2, and intro- 
duce a wave function by separating the uniform motion 
of the bound particle as a whole, from the internal mo- 
tion: < E > p(.t) = e lPX x(xi, x 2 ). Finally we obtain: 

$ P (x) =-2 v / ^S'(x)7 5 sin(Pa;/2) 

+ -^e 7E e l92/3(x) e- je75 7 5 cos( J Px/2) . (37) 

2V7T 

This amplitude wave function entirely agrees with that 
found in (|24j) from the factorizable part of the bound 
state residue. 



VI. SUMMARY 

In conclusion we would like to recapitulate our main 
results. We have found analytically the polar contribu- 
tion to the qq Green's function. In the topologically triv- 
ial case, the residue turns out to be factorizable, which is 
commonly taken for granted. Out of this we obtained the 
formula for the Bethe-Salpeter wave function. The situa- 
tion turned out to be different, if one considers one or two 
instanton background. We have obtained expression for 
the residue in this case and found, that the BS function 
has been modified by the 0-vacuum structure. The full 
polar term may now be written as a sum of several fac- 
torizable terms. One sees that, apart from the formation 
of the bound state (Schwinger boson), there appears also 
another way for a pair qq, not bearing meson quantum 
numbers, to travel over the space. They may jump from 
one instanton to the other. This phenomenon is possible 
in there are at least two instantons. 

Moreover if one takes into account also these compo- 
nents, that are usually neglected for P 2 — > as being 
small in comparison with that revealing a pole, one finds 
another type of singularity — a branch point connected 
with the appearance of logarithmic and dilogarithmic 
functions. They arise because of long-range correlations 
in the condensed vacuum. 

All these instanton contributions vanish, if we consider 
the 5-matrix elements, i.e. only transition amplitudes for 
the 'physical' asymptotic states. These, in the Schwinger 
Model, are Fock states of one or more bosons. The four- 
point (quark) function is then replaced with the current- 
current function, for which it is shown, that simple pole 
structure as well as factorization of the residue reappear. 
The analytical structure of the Green's function may then 
be more intricate than that of the S-matrix, especially for 
topologically nontrivial theory. 
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